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I. INTRODUCTION 



Specific aspects of hadronic structure are described by several different phenomenological functions. Form factors, 
usual parton distribution functions (PDFs) and distribution amplitudes are the so-called "old" phenomenological 
functions since they have been around for a long time. On the other hand, the concept of generalized parton 
distributions (GPDs) Q, IM IM |^ (for recent reviews, see 0, 01) is new. These "new" phenomenological functions 
are hybrids of the "old" ones, and therefore they provide a unified and more detailed description of the hadronic 
structure. 

In recent years, significant effort was made to access GPDs through the measurement of hard exclusive electro- 
production processes. The simplest process in this respect is the deeply virtual Compton scattering (DVCS) process, 
i.e. an electron scatters off a nucleon producing a photon in the final state. Using the neutrino beam instead we 
extend the DVCS process into the weak sector, where one expects to be sensitive to a different fiavor decomposition 
of GPDs, and further, due to the presence of the axial part of the V-A interaction also sensitive to a different 
set of GPDs. In addition, the weak DVCS process allows us to study fiavor non-diagonal GPDs, e.g. in the 
neutron-to-proton transition. 

Detailed study of weak deeply virtual Compton scattering with electron and neutrino beams will be presented in a 
forthcoming paper. 



II. PHENOMENOLOGICAL FUNCTIONS 



Among phenomenological functions we discuss form factors, usual parton distribution functions and generalized 
parton distributions. 



A. Form Factors 



Form factors are defined through matrix elements of electromagnetic and weak currents between the hadronic 
states. In particular, the matrix element of the vector current between the nucleon states N (pi, si) and N {p2, S2) is 
parametrized in terms of the nucleon electromagnetic form factors, i.e. Dirac and PauH form factors. 



2M 



(1) 



where r — pi — p2 is the overall momentum transfer, the invariant t ^ r"^ and M denotes the nucleon mass. In the 
case of the axial vector current one has the axial and pseudoscalar form factors. 



{N{p2,S2)\J'XiO)\N{p^,S^)) - U{P2,S2) 



u{pi,si). 



(2) 



Both currents in Eqs. Q and Q are given by the sum of their fiavor components, Jy (0) = J2f Qfi^f {^)l^'^f (0) 
and (0) — J2f Qaj^j {Q)l^lb'4'f (0), where Qj and Qaj are the electric (in units of |e|) and axial charges of the 
quark of fiavor /, respectively. A similar decomposition holds for the form factors, Fi^2 (t) — J2f QfFi,2f (t) and 
9A,p (t) = J2f QAfdA.Pf (t). Their limiting values at i = are known, e.g. Dirac and Pauli form factors give the 
total electric charge of the nucleon and its anomalous magnetic moment. 



In particular, the nucleon electromagnetic form factors can be measured through elastic electron-nucleon scattering. 



The process is shown in one photon-exchange approximation in Figure left. 



(3) 



3 



B. Usual Parton Distribution Functions 

PDFs are defined through forward matrix elements of quark/gluon fields separated by light-like distances. For the 
unpolarized case one has 



{N (pi, Vi/ i-z/2)r^f (z/2) |iV (pi, 



u(pi,si)7^u(pi,si) 



X I dx 

10 



-ix{piz) 



(4) 



and for the polarized one 



{N (pi, si)| -ipf {-z/2) 7^75-0/ \N (pi, si))^2^o = u (pi, si) 7''75W (pi, si) 



X dx 
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e-*^(Pi^)Aq^ (x) + e"(Pi^)Ag^ (x)] . (5) 



To make connection with GPDs, which are usually discussed in the region —l<x< 1, it is convenient to introduce 
new distribution functions, 



Qf (x) = 



Qf (x) X > 

-qj {—x) X <Q 



and 



(6) 



^If - I Ag^ i-x) X < 
and alternatively, write the integrals over x in Eqs. Q and JHl as 



(7) 



dx 



-ix(pxz] 



If (x) 



Jx(piz] 



dx 



-ix{piz) 



Aqf (x)+e"(Pi")Ag^ (, 



dx e-"(Pi^)g/ (a;) , 
dx e-'''^P'''> Aqf (x) . 



(8) 



Furthermore, one observes that the definition of PDFs has the form of the plane wave decomposition. Thus it allows 
us to give the momentum space interpretation: qf (x) /qf (x) is the probability to find the quark/antiquark of fiavor 
/ carrying the momentum xpi inside a fast-moving nucleon having the momentum pi . 

PDFs have been intensively studied in hard inclusive processes for the last three decades. The classic example in 
this respect is the deeply inelastic scattering (DIS) process, i.e. inclusive scattering of high energy leptons on the 
nucleon, 



e'N — > e-X, (9) 

shown in Figure^] right. It played a key role in revealing the quark structure of the nucleon. The structure functions, 
accessed in the DIS process, are directly expressed in terms of PDFs. Through the optical theorem, its cross section 
is given by the imaginary part of the forward virtual Compton scattering amplitude (see Figure ^ . The summation 
over X refiects the inclusive nature of the nucleon structure description by PDFs. 

In the Bjorken regime, where the space-like momentum transfer is sufficiently large together with large total center- 
of-mass energy of the photon- nucleon system, —qf and (pi + gi) oo, while the ratio xb = — 9i/2 (pi ■ qi) is finite, 
perturbative QCD factorization works. In other words, the forward virtual Compton scattering amplitude factorizes 
into perturbatively calculable hard scattering process at the level of quarks and gluons, and process independent 
matrix elements which contain the soft non-perturbative information about the nucleon structure represented by 
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the blob (see Figure We recall that these forward matrix elements consist of quark and gluon operators, whose 
fields are separated by a Hght-Hke distance. They are described and parametrized in terms of PDFs. Schematically, 
QCD factorization allows us to write the amplitude in the form of the so-called handbag diagrams. Moreover, the 
leading contribution in the lowest order in the strong coupling constant as is given by two (s- and tt-channel) handbag 
diagrams in which the hard propagator is convoluted with the PDFs. Taking the imaginary part of the forward virtual 



Compton scattering amplitude generates the delta function, 3 

the DIS process one measures the PDF qf {x) at two points, x 
and X — —xb for that of antiquarks. 



{xpi ±qiy + it 
= ±XB, with X 



~ 5 (x =F xb) /2 {pi ■ qi). Hence in 
XB corresponding to the quark PDF 



C. Generalized Parton Distributions 



A more recent attempt to use perturbative QCD to extract new information about the hadronic structure is the 
study of hard exclusive electro-production processes, in particular the DVCS process. This is a much more difficult 
task due to the small cross sections. However, high energy and high luminosity electron accelerators combined with 
large acceptance spectrometers give a unique opportunity to perform precision studies of such reactions. The DVCS 
process can be accessed through the reaction 



e^N — > e^N-f. (10) 

It turns out that factorization into short and long distance dynamics is more general. Having large space-Hke 
virtuality of the initial photon while the final state photon is on shell is sufficient for QCD factorization to work. In 
contrast to the DIS process, the outgoing photon is real, and henceforth, the overall momentum transfer is not equal 
to zero. In the leading handbag approximation, the so-called non-forward virtual Compton scattering amplitude is 
dominated by two diagrams (see Figure 21 • The lower blob now contains the non-forward matrix elements of the 
same quark and gluon operators as in the forward case. They are parametrized in terms of GPDs. 



It is convenient to introduce the average of the nucleon momenta, p — (pi +P2) /2, and treat the initial and 
final hadron in a symmetric way (see Figure EJ. In this scheme, at the leading twist-2 level, the nucleon structure 
information can be parametrized in terms of two unpolarized and two polarized GPDs denoted by H, E and iJ, 
respectively. They are functions of three variables (a;,^,t), and further they are defined for each quark fiavor /. In 
addition to the usual light-cone momentum fraction x, GPDs also depend upon another scaling variable, the skewness 
parameter ^ = r||/2p||, specifying the longitudinal momentum asymmetry, and upon the invariant t. The variables x 
and ^ solely characterize the longitudinal momenta of the partons involved, however, the i-dependence of GPDs is 
related to their transverse momenta. Thus one can simultaneously access the longitudinal momentum and transverse 
position of the parton in the infinite momentum frame Furthermore, by removing the quark with the light-cone 
momentum fraction a; -I- ^ and replacing it with the quark of the momentum fraction x — ^, one can say that GPDs 
measure the coherence between two different quark momentum states of the nucleon, i.e. the quark momentum 
correlations in the nucleon, whereas usual PDFs yield only the probability that a quark carries a fraction x of the 
nucleon momentum. 

Since — 1 < a; < 1 in FigureEl the momentum fractions x±£^ of the active quarks can be either positive or negative. 
Positive and negative momentum fractions corresponds to quarks and antiquarks, respectively. Therefore GPDs 
have three distinct regions: when £, <x <1 (— l<a;< — ^) both partons represent quarks (antiquarks), whereas for 
—i<x<£, one parton represents a quark, and the other parton an antiquark. In the first two regions GPDs are just 
the generalizations of the usual PDFs, however, in the third region they behave Hke a meson distribution amplitude. 
Hence, in the region — ^ < a; < ^, they contain new information about the nucleon structure since this region is not 
present in DIS. 

GPDs have interesting properties Hnking them to usual PDFs and form factors. In the forward limit, pi = p2 and 
r, t = 0, the GPDs H and H coincide with the quark density distribution qj {x) and the quark helicity distribution 

Ag/ {x) given by Eqs. © and obtained from the DIS process. One writes the so-called reduction formulas for the 

functions H and H , 



i7/(x,0,0) = |_J 
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and 

while the functions E and E have no connections to PDFs. They are always accompanied with the momentum transfer 
r, and therefore invisible in inclusive measurements. In the local limit, z — 0, GPDs reduce to the form factors. In 
other words, the first moments of GPDs are equal to the nucleon elastic form factors. Namely, 

dx Hf{x,^,t)^Fif{t) , / dx Ef{x,^,t)^F2f{t), 
-1 J-i 

dx Hf{x,tt)=gAf{t) , / dxEf{x,^,t)=gpf{t). (13) 
'-1 J-i 

We call these relations the sum rules. They are model and ^-independent. 

GPDs are also relevant for the nucleon spin structure. In particular, the second moment of the unpolarized GPDs 
aX t — gives the quark angular momentum. 



1 

= 2^ dxx[Hf{x,^,t^O) + Ef{x,tt^O)]. 



(14) 



The above equation is independent of The quark angular momentum, on the other hand, decomposes into the 
quark intrinsic spin and the quark orbital angular momentum. 



J, = -^^ + L,, (15) 

where AS is measured through the polarized DIS process. Substituting Eq. lfn)l into Eq. lfT5|l one can determine Lq. 
Moreover, since the nucleon spin comes from quarks and gluons, 1/2 = Jq + Jg, one can extract the gluon contribution 
Jg to the nucleon spin. Hence by measuring GPDs one obtains information about the angular momentum distributions 
of quarks and gluons in the hadron. 



III. WEAK DVCS AMPLITUDE 



In the most general case, the virtual Compton scattering amplitude is given by a Fourier transform of the correlation 
function of two electroweak currents. In particular, for the weak DVCS process we have 

T^" = ij d^ze-^'>^N{p-r/2,S2)\T{Ji^{zl2)rE,j{~zl2)}\N{p + Tl2,s^)), (16) 

where corresponds to either the weak neutral current J^jy (-^Z^) or to the weak charged current {z/2), 

i.e. to the exchange of the weak boson or W^, respectively. One of the methods to study the behavior of Eq. 
ifTfill in the generalized Bjorken region is to use the light-cone expansion for the time-ordered product of two currents 
T {j!^ {z/2) J'^j^j {—z/2)} in the coordinate representation. The expansion is performed in terms of QCD string 
operators. Its leading order contribution is shown in Figure El The hard part of both handbag diagrams starts at the 
zeroth order in as with the purely tree level diagrams in which the weak virtual boson and real photon interact with 
quarks being treated as massless fermions. Since the weak current couples to the quark current through two types of 
vertices, qqZ^ and qqW^, the quark fields at coordinates ±z/2 can carry either the same or different flavor quantum 
numbers. We treat these two cases separately. 
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A. Weak Neutral Current 

We expand the time-ordered product of the weak neutral and electromagnetic current in Eq. Ijlfijl , 

iT{J^^{z/2)rEM{-zm} = -^5]Q/[vi/(-^/2)7V7''(4-75C^)V'/(^/2) 

-i>f (z/2) 7^ (4 - 75C^) Yl'^f (-^/2)] , (17) 

f f 

where Cy and c\ are the weak vector and axial vector charges, respectively. We express the original bilocal quark 
operators with three Lorentz indices in Eq. itTTIl in terms of the string operators with only one Lorentz index, 

^TUi^r, (^/2) Jem (-^/2)} = E Qf {'^v'""'' (-^Z^) ^n^f (^/2) " ^ -^)] 

TT Z ^ 

+cite^'"'^ [4,f (-z/2) 7,75V/ {z/2) + (z — > -z)] 
-c^s^""'' [i^f (-z/2)7,75V/ {zm - (z ^ -z)] 

-c^ze'''"-" (-Z/2) 7„ V/ (z/2) + (z ^ -z)] } . (18) 

We have two types of string operators in Eq. ifTsji . The (axial) vector string operators come (with) without 75. 
They can be accompanied with tensors, s'^'""' and e'^'""', which are symmetric or antisymmetric in indices ^, v. 
Furthermore, in addition to the standard electromagnetic DVCS process, we end up with two more terms. 

In the next step we isolate the twist-2 part of the string operators in Eq. (jl8ll . and sandwich it between the 
initial and final nucleon state. At this point we introduce the relevant non-perturbative functions (GPDs) in order to 
parametrize the non-forward nucleon matrix elements of the vector and axial vector string operators on the light-cone. 
Namely, 

(A^(p2,S2)|^/(-z/2) /V/(z/2)±(z^-z)|iV(pi,si))^,^o = u{p2.S2) Mpi,si) J ^dxe-"^P'^Hfix,^,t) 

+U(P2,S2) — m(P1,Si) 



{N ip2,S2)\^f {~z/2) /j5^f {z/2)T (z ~z)\N {pi,si))^,^„ = u(p2, S2) ^75u(pi, si) / dx e-'^^^'^ Hf {x,^,t) 



-w(P2,S2)-^^75m(pi,si) 

X j dxe~"^P'')Ef{x,^,t). (19) 

The "plus" distributions enter the electromagnetic DVCS process, and they correspond to the sum of quark and 
antiquark contributions, i.e. to the sum of contributions from valence quarks and twice the sea quarks. On the other 
hand, the weak version of the process gives access to the "minus" GPDs, which correspond to the difference of quark 
and antiquark contributions. This difference is equal to the valence quark contribution. 

B. Weak Charged Current 

The expansion of the time-ordered product of two currents in the weak charged sector reads 

^T{J^^iz/2)J'^,,i-z/2)} = --^^[Q^,vV(-^/2)7V7'' (1-75)^-/(^/2) 

/,/' 



-Qf{Pf (z/2) r (1 - 75) 7^7'' V/' i-z/2)] . (20) 
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Here the sum over quark flavors is subject to an extra condition, Q f — Q f = 1 or —1, due the fact that the weak 
virtual boson carries an electric charge. Hence the initial and final nucleon are not the same particles anymore. 
Now the vector and axial vector string operators obtained from Eq. lf2(Hl are accompanied by different electric charges 
and quark flavors. For that reason the non-forward nucleon matrix elements, 

{N ip2, S2)\ (z) \N (pi, Si)) - {N {p2, S2)\ (-^/2) (^/2) ±{f'^f,z^ -z) \N (pi, si)) , 
{N{p2,S2)\di'Uz)\N{p,,si)) = (7V(p2,S2)|Vi/'(-2/2) /75^/(z/2)T(/'^/,^^-2)|iV(pi,si)), (21) 

involve different flavor combinations. They are parametrized in terms of GPDs, which are non-diagonal in quark 
flavor, 

{Nip2,S2)\df'^ iz)\N{p,,s,))^,^, 



{Nip2,.S2)\d{Uz)\N{p,,s^))^.^, 



IV. WEAK DVCS PROCESSES 

Before we examine speciflc processes we introduce a simple model and discuss the kinematics which is common to 
all DVCS-like reactions. 

Our simple model has three properties. First we assume that the sea quark contribution is negligible, and therefore 
the "plus" GPDs are equal to the "minus" GPDs with the quark flavor f = u,d. Secondly we take a factorized ansatz 
of the i-dependence from the other two scaling variables x and ^ for all distributions. The i-dependence of GPDs is 
characterized by the corresponding form factors given by Eq. Thirdly we neglect the ^-dependence in all GPDs 

except in the E distribution. The parametrization of GPDs is taken from Refs. P. fl(iL ITlL fl^ . Namely, for the H 
and E distributions one has 

ffr' C, t) = qT' [x] and H^' {x, t) = g^' {x) E,, (t) , 

ET' {x, 1 1) = qT' [x) and ET' (x, 1 1) - qT' (t) , (23) 

where the unpolarized valence quark distributions in the proton are given by 

(x) = 1.89a;-°-* (1 - xf"^ (1 + 6x) and (x) = 0.54x-°-'^ (1 - xf^ (1 + 82;) . (24) 

The u- and d-quark form factors can be extracted from the proton and neutron form factors, which can further be 
related to the Sachs electric and magnetic form factors. For the polarized quark GPD H one has 

i?r' t) = ix) f 1 + ^) and H^' (^, C, t) = Aq^' (x) f 1 + ^) , (25) 

with the mass m 4 = 1.03 GeV. The polarized valence quark distributions can be expressed in terms of the unpolarized 
ones through 

Agr' = cos9d (qT' - IqT') and AqT' = COS 0c ( -W') , (26) 



+U{P2,S2)^-^^^^U 



■{P2,S2) Mpi,si) I dxe-'^'^P'^Hf,f{x,^,t) 



-U{p2 



,S2)^75^^(pi,si) l\ dx e-^-^^^^Ef,^ {x,^,t) . 



(22) 
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where cosBd = [l + Ho (l — x^) /v^] ^^'^ Hq = 0.06. Finally, for E we accept the pion pole dominated ansatz 

^r'(x,e,i). (27) 



(t) = AgA [t = 0) M2 

where denotes the pion mass and qa [t — Q) — 1.267. 

In general, a neutrino induced DVCS process on a nucleon is given by the reaction 

/i(fc)+iVi(pi) l2{k')+N2{p2)+l{q2), (29) 

where a neutrino li scatters from a nucleon A^i to a final state I2, nucleon N2 and real photon 7. Schematically, 
the reaction is presented by three diagrams (see Figured- The first diagram is the so-called DVCS diagram which 
corresponds to the emission of the real photon from the nucleon blob. In our approximation it is calculated from two 
handbag diagrams (see Figure ■ In the other two (Bethe-Heitler) diagrams of Figure the real photon is emitted 
from a lepton leg. 

The differential cross section in the target rest frame, in which the weak virtual boson four-momentum qi has no 
transverse components, assumes the form 

d'^ ^ J 1 1 + xb{M/lo) 2 

dxBdQjdtdip 64:3 (2TTf Muxsiy + 2xB{M/uj)f 

where T represents the invariant matrix element. The invariants in Eq. (13(1)1 are given by s = {k + pi)'^ , Q\ = 
—ql, and y = (pi ■ qi) / (j>i ■ k). Moreover, lj denotes the energy of the incoming neutrino beam, and ip the angle 
between the lepton and nucleon scattering planes. One finds the kinematically allowed region for the reaction (see 
FigurelHll under the following constraints: fixed neutrino beam energy at w = 5.75 GeV, the invariant mass squared of 
the weak virtual boson-nucleon system, s = {pi + qif > 4 GeV^, and the virtuality of the weak boson, Ql>2 GeV^. 
As an example of the particular kinematics, we further choose Ql = 2.5 GeV^, xb = 0.4 and set ip = 0. One plots 
i as a function of the angle Bbj between the incoming weak virtual boson B and the outgoing real photon 7 in the 
target rest frame (see Figure inj. 

We study two different examples of the weak DVCS processes. In the weak neutral sector we consider neutrino 
scattering off an unpolarized proton target through the exchange of Here one only measures the Compton 
contribution (see Figure since the photon can not be emitted by the neutrino. Next we consider neutrino-neutron 
scattering through the exchange of W'^ with a proton in the final state. This particular process, however, has both 
contributions (see Figure El • The Bethe-Heitler background is calculated from one diagram (see diagram (b) in 
Figure since only the outgoing muon can emit the real photon. In contrast to the standard electromagnetic DVCS 
process on an unpolarized proton target (see Figure [l2)l . where the Bethe-Heitler cross section is well above the 
Compton one, the situation in the weak charged sector is the other way round. 

V. CONCLUSIONS 

Generalized parton distributions contain the most complete and unified description of the internal quark-gluon 
structure of hadrons. Form factors, usual parton distribution functions and distribution amplitudes, on the other 
hand, can be treated just as particular or limiting cases of generalized parton distributions. Furthermore, the 
formalism of generalized parton distributions provides nontrivial relations between exclusive and inclusive processes. 



ET' (x, t) = (t) 0^ j and Er (x, t) 

The function F^r (t) and the pion distribution amplitude (pT^ are taken in the form 



1.7 



(TO2+|t|)/GeV2 (l + |t|/2GcV2)' 



and (j)TT {u) = 6m (1 — m) , 



(28) 
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and also between different exclusive processes. 

We have extended the deeply virtual Compton scattering process into the weak sector by using the neutrino beam. 
We have argued that the weak deeply virtual Compton scattering process gives an additional information about the 
hadronic structure. It is expected in the near future that neutrino scattering off a nucleon will be studied at high 
intensity neutrino beam facilities. 
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VI. FIGURES AND PLOTS 




FIG. 1: Elastic electron-nucleon scattering in one-photon exchange approximation (left) and deeply inelastic scattering (right). 




FIG. 2: Optical theorem. 




FIG. 6: Handbag diagrams (s- and u-channel) in the weak DVCS amplitude. 
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a b c 



FIG. 7: DVCS (a) and Bethe-Heitler (b and c) diagrams contributing to the lepton-production of a real photon. 
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FIG. 8: Kinematically allowed region. 




FIG. 9: Momentum transfer squared plotted as a function of the angle between the incoming weak virtual boson and outgoing 
real photon in the target rest frame. 
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FIG. 10: Weak neutral DVCS differential cross section on an unpolarized proton target plotted as a function of the angle 
between the incoming weak virtual boson and outgoing real photon in the target rest frame. 
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FIG. 11: Compton (dashed line) and Bethe-Heitler (solid line) contribution to the weak charged DVCS differential cross section 
on an unpolarized neutron target plotted as a function of the angle between the incoming weak virtual boson and outgoing real 
photon in the target rest frame. 




FIG. 12: Compton (dashed line) and Bcthe-Heitler (solid line) contribution to the electromagnetic DVCS differential cross 
section on an unpolarized proton target plotted as a function of the angle between the incoming weak virtual boson and 
outgoing real photon in the target rest frame. 



